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We study the application of AdS/CFT duality to longitudinal boost invariant Bjorken expan- 
sion of QCD matter produced in ultrarelativistic heavy ion collisions. As the exact (1 + 4)- 
dimensional bulk solutions for the (1 + 3)-dimensional boundary theory are not known, we 
investigate in detail the (1 + l)-dimensional boundary theory, where the bulk is AdSs grav- 
ity. We find an exact bulk solution, show that this solution describes part of the spinless 
Bahados-Teitelboim-Zanelli (BTZ) black hole with the angular dimension unwrapped, and 
use the thermodynamics of the BTZ hole to recover the time-dependent temperature and 
entropy density on the boundary. After separating from the holographic energy-momentum 
tensor a vacuum contribution, given by the extremal black hole limit in the bulk, we find that 
the boundary fluid is an ideal gas in local thermal equilibrium. Including angular momentum 
in the bulk gives a boundary flow that is boost invariant but has a nonzero longitudinal 
velocity with respect to the Bjorken expansion. 
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1 Introduction 



Collisions of large nuclei at very high energies can be modelled by taking the transverse size 
and collision energy to be effectively infinite, so that the dynamics is invariant under boosts 
in the longitudinal direction and translations in the two transverse directions. When the 
Minkowski metric is written as ds^ = —dt^ + dx^ + dx2 + dxg, where x is the longitudinal 
direction, the outcome of the collision takes place in the wedge t > the natural coordinates 
are r = \Jt'^ — x'^ and ?? = (1/2) log[(t+3;)/(t— x)], the metric reads ds^ = —dT^+T^dri'^+dx2 + 
dx"^, and the hydrodynamic variables are independent of rj, X2 and X3. Denoting the shear 
and bulk viscosities by respectively fj and (, the hydrodynamic equations read V ^T^^ = 0, 
T^^^ = (e + + pg^^^^ + AT'^'^, AT^'' = + Oid'"'' - ui^u"), e{T) = 3p{T) = SaT^, 

fj = p'(T)/(47r) = aT^/n, a = ir'^N^ /8, and C = 0- For the longitudinal similarity flow, 

x) = — = tanh 0(r, ?]), @{T,rj)=r], = — , (1) 

t T 

these equations can be directly integrated to give 

1 \ rr,v/' 1 



T{t) = r, + 



GlTTi J \T J GlTT 



where Tj is normalised by T(Tj) = Tj and all the constants refer to the by now standard 
predictions for = 4 supersymmetric gauge theory [HIS]- Determining these initial values 
for colliding nuclei is an important problem in QCD dynamics [3]; very approximately, they 
are 0] Ti,Ti = 0.5 GeV, 0.2 fm/c at ^/s = 200 GeV (Relativistic Heavy Ion Collider energies) 
and Ti,Ti = l GeV, 0.1 fm/c at = 5500 GeV (CERN LHC energies). 

A study of expanding systems in the gauge theory/gravity duality picture was initiated in 
O El E] and continued in several further papers P (S HOl EJ HSJ [ill [Ml . In [7j the starting 
point was to write a candidate five-dimensional gravity dual of the above collision process in 
the form 

£2 



ds = ^ 



-a{T, z)dT^ + T^h{T, z)dr]^ + c(r, z){dx2 + dx^) + dz^ (3) 



and then study what constraints five-dimensional Einstein's equations give for the functions 
a, b and c and to the holographic energy-momentum tensor computed from them. In par- 
ticular, one may expect to measure the last term in ([2|) and thus independently check the 
standard prediction rj/s = l/(47r). However, since the holographic energy momentum tensor 
determines only e ~ T^, not T directly, one is not able to measure the last term in ([2]) inde- 
pendently, and interference terms will appear, even if one relates energy density directly to 
the viscosity |12j . Also, no exact solution for a, b and c is known and it is difficult to judge 
the validity of the several interesting results obtained by considering the large-r behavior of 
the solutions. One may further ask whether the time dependence in ([3]) could be removed by 
a coordinate transformation, as the case is for metrics admitting the larger isometry group 
-B3 [15] . We shall in the following simplify the problem by neglecting entirely the transverse 
dimensions X2 and X3. In the this case Einstein's equations imply that the bulk geometry is 
locally AdSs, and the bulk side becomes exactly solvable. 
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In case of AdSs the appropriate 1-brane solution [16] in lOd combines Qi Dl-branes in 
the X5 direction with Q5 D5-branes in the x^, ...,x^ directions. The x^, ...,x^ directions are 
compactified on a 4-torus T4 with V4 ~ a'^ and the x^ is taken to be of length L ^ \/^. 
When the x^, are written in spherical coordinates with ds'^ = dr"^ + r'^d^\, the 0, 5,r 

coordinates in the r — > limit give us AdSa and the whole structure is AdSsxSa x T4. The 
boundary theory dual to string theory in this background is a 2d conformal field theory in 
x^,x^ with 4Qi(55 bosons and an equal number of fermions. We shall neglect the dilaton 
and the 3-form field strength and determine the metric by solving AdSs gravity equations. 
Several different coordinate systems are studied. The holographic energy momentum tensor 
is determined. We observe that the energy density computed in this way is exactly that of 
an ideal gas of ^QiQ^ bosons and fermions, the factor 3/4 observed for the usual AdSs case 
is 1 now. 

The differences in the application of AdS/CFT duality to AdSsxSs and AdSsxSs x are 
manifest in the relation between string theory and supergravity background parameters. For 
AdSsxSs we have 

C' = A7:g,N,a'^ = g\^N,a'\ ^ = (4) 

G5 vr 

where the string coupling gg is constant since the dilaton is a constant. Thus the AdSs radius 
and the string tension are simply related via the coupling constant gYM of the boundary 
theory, M = A super- Yang-Mills (SYM). For AdSsxSs x 

= gl^^^QiQ,a'\ ^ = 4QiQ5. (5) 

Thus the relation between the string tension and AdSa depends on the compactification 
volume V4, which is not an experimental number. However, for both cases the dimensionless 
relation between C and the Newton constant is very simple, ~ number of degrees of freedom. 

While the integration of hydro dynamical equations V^T'^'^ = is trivial for given initial 
conditions, the real problem is in the determination of initial conditions. For heavy ion colli- 
sions these will depend on the atomic number A of the colliding nuclei and the collision energy 
[31 S] . AdS / CFT even in the well controlled case of = 4 SYM permits the determination 
of vacuum expectation values of operators with fields integrated over —00 < t < 00, but the 
application to path integrals with fields starting at some ti seems to be an open question. 



2 The bulk solution 

The local properties of solutions to AdSs gravity equations are well known [17^ [T8] , but since 
we aim at a solution with a specific structure, the derivation gives some insight. The AdSa 
action with the cosmological constant A = 1 and the 3d Einstein equations are 

1 . 2\ 



d-'x^{R+-. , (6) 



Rab — -^Rdab — -^9ab — 0, (7) 
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which further imply 
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As the (2 + l)-dimensional analogy of ([3|), we adopt the ansat^ 

£^ r 



(8) 
(9) 

(10) 



where z > and — oo < 77 < oo. The nontrivial components of Einstein's equations then 
yield the following set of equations: 

dztt — zdia = 0, 
d^b - zdlb = 0, 
dza drb = a dzdrb, 



TT : 
rjri : 
TZ : 



zz : —a^dzb + a^dza{—b + zdzb) + zdradrb — za d^b = 

We look for solutions that have at z — > the asymptotically AdS form [19] 

£2 



(11) 
(12) 

(13) 

(14) 



where the two-dimensional metric g^y has the small z expansion 
gAr.z) = 5?J(t) +5?J(r)z2 +5(t)(r)z^ + 
and the conformal boundary metric gl^J is the Milne metric |20j 



gj^Jdx^'dx'' 



-dr"^ + T'^drf, 



(15) 
(16) 
(17) 



with < T < 00. Equations (jlip and (112p integrate immediately to 

a(T,z) = ai(r)z2 + a2(T), (18) 
5(t,z) = 6i(t)z2 + 62(t). (19) 

Matching to the boundary metric (jl7p gives a2(''") ~ 1 ^"^^ ^K''") = Equations (|13p and 
(jl4p then reduce to a form from which ai and 61 can be found by elementary integration. 
The general solution can be written as 



ds' 



z2 



1 



(M - l)z^ 
4^2 



dT'+\l+ 



(M - l)z2 
4^2 



T'^dif 



dz^ 



(20) 



where the constant M may take any real value. 

We shall from now on assume M to be positive. It will be explained at the end of section 
[3] why this is the only case in which we can recover interesting thermodynamics. 



^One can also start from an ansatz with 77-dependence, i.e. a{T,z) a{T,ri,z) and 6(r, 2) b{T,Tj,z), but 
if one wants the diag( — 1, r^) boundary metric one can show that all rj dependence dies away. The absence of 
cross-terms in the ansatz is crucial for this result. 
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3 Relation to the AdSs black hole 



The metric (|20p is explicitly time dependent, and it has a coordinate singularity at = 
\M — l\z^ / A:. As a solution to Einstein's equations ([7]), the metric must however be locally 
AdSs |17l I18j . We shall now show that the metric covers part of the spinless Bahados- 
Teitelboim-Zanelli (BTZ) black hole whose angular dimension has been unwrapped, and the 
coordinate singularity resides in the white hole region of this spacetime. 
Introducing the coordinates (f7, V) by 



2t 



U 



V 



2t+ (/M + l) ^ 
'2r - f ^/M + i)z 



2t+ I^/M - 1) z 




(21) 



the metric (1201) takes the form 



AdUdV ,^ /I - UVV , 



(22) 



Suppose for the moment that rj were periodic with period 27r. The metric (|22|) is then the 
spinless nonextremal BTZ black hole, and the global Kruskal-type null coordinates (U, V) 
have the range — 1 < UV < 1 [n\ I18j . The left and right infinities are at UV — > —1, the 
Killing horizon of the Killing vector Vdy — Udjj is at UV = 0, and the black and white hole 
orbifold singularities are at UV 1. The conformal diagram is shown in Figure [TJ In each 
of the quadrants in which UV ^ 0, the coordinate transformation 



l-UV 
1 + UV 



2VM 



In 



(23) 



brings the metric to the Schwarzschild-like form 



ds' 



M 



dt^ + 



dr'^ 



r'^jC? - M 



+ r^dr^, 



(24) 



in which dt = {yMj C^jiydy — Udjj), and the Killing horizon of dt is at r = LyM. The 
Arnowitt-Deser-Misner (ADM) energy equals Mj (8G3), where the factor 1/(8^3) comes from 
our normalisation of the Einstein action . 

In our spacetime r\ is not periodic, and the singularities of ([22|) at UV 1 are just 
singularities of the coordinate system {U,V,r]) on the AdSs hyperboloid [171 HI]- What is 
important for us is that the observations about the Killing vector dt remain valid even when 
r] is not periodic. We shall continue to describe the metric (j22p in the black hole terminology 
even for nonperiodic r], trusting that no ambiguity will ensue. 
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r = 



Figure 1: The conformal diagram of tlie spinless nonextremal BTZ black-and-wliite hole ([22]) [171 [18]. 
The coordinate rj is suppressed, the coordinates of the diagram are p = arctan U and q = arctan V 
with |p + ql < tt/2 and |p — f/l < 7r/2, and the p-axis (respectively q-axis) is tilted 45 degrees to the 
left (right) from the vertical. The Killing horizon of the Killing vector dt = {sfM / L){Vdv — Udu) is 
at pq = 0, where r = r+ = C^/M . The coordinate grid shows the curves of constant r (solid curves) 
and the curves of constant z (dotted curves) in the region < z < wr, where v is given by (|25p and 
Af > 1 is assumed. The dashed line is the coordinate singularity z = vt, r = M — 1. On a given 
curve of constant r, z takes the values < z < vt, increasing from right to left, and on a given curve 
of constant z, t takes the values z/v < t < oo, increasing bottom to top. 



The region of the BTZ hole ([22]) covered by ([20]) depends on whether < M < 1 or 
1 < M < oo. Suppose first that 1 < M < oo. The metric (I20p is then regular for z < vt and 
VT < z, where 

2 , ^ 

v= , 25 

Examination of ([2ip shows that the region z < vt covers in ([22p the region 



U<0, V>{ ^ - , (26) 




and the coordinate singularity at z = vt is on the spacelike surface of constant UV where the 
latter inequality in ([26p becomes saturated. The curves of constant r, the curves of constant 
z and the coordinate singularity are shown in Figure [TJ If we take r to increase to the future, 
so that U and V increase to the future, the region z < vt thus comprises one exterior region, 
where C\fM < r < oo, and the part £\/ M — 1 < r < L\fM of the white hole interior 
region. The coordinate singularity is on the spacelike surface r = LyjM — 1 in the white hole 
interior. The parameter v has an interpretation as the coordinate velocity of the coordinate 
singularity, but this coordinate singularity is neither an event horizon, apparent horizon nor 
a dynamical horizon [211 [22], and would not be even if r] were made periodic to satisfy certain 
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technical conditions in the definitions of these horizons. The coordinate singularity arises 
just because the surfaces of constant r become parallel to the surfaces of constant z/r. 

Similar considerations hold for the region vt < z, and the region covered in (|22p is obtained 
by interchanging U and V in ()26p . A qualitative picture of the curves of constant r and the 
curves of constant z is obtained by a left-right interchange in Figure [TJ As we are interested 
in the conformal boundary where z ^ with fixed r and rj, the region of interest for us is 
however z < vt. 

Suppose then that < M < 1. The coordinate singularity is now at r = — M z, 
and in the conformal diagram of Figure [1] it is at the white hole singularity. The region 
— M z < T, which reaches the conformal boundary at z ^ 0, thus covers in (j22p one 
exterior and all of the white hole interior. 

For use in section U we record here some thermodynamical properties of the BTZ metric. 
For quantum fields in AdSs, an AdS-invariant vacuum state induces in the metric (j24p a state 
that is thermal with respect to the Killing vector dt, in the BTZ temperature [17\ [T8] 



The ADM energy AE and the Bekenstein-Hawking entropy ASbtz for an interval Arj are 

A^; M ASbtz Vm 



At] leTrGg' CAr] AG3 



(28) 



These formulas originally arose in the context in which rj has period 2tt and the bulk spacetime 
is the BTZ black hole, but they remain valid also for nonperiodic t] if one maintains that all 
thermodynamical quantities should be defined so that they are invariant under translations 
in 77. It is this assumption of translational invariance that forces one to regard dt in (j24p as the 
physically relevant time translation Killing vector, as the only Killing vectors that commute 
with translations in rj are 5^ and dt- For the BTZ hole the translational invariance comes 
from the periodicity of rj, while for us it is motivated by the aim to model a boost-invariant 
flow on the conformal boundary. 

For use in section [H we also recall that the coordinate transformation 

(where z is not the same as that in (|2Up ) brings the metric (|24p to the asymptotically AdS 
standard form of (1151) and (I16D. 



£2 

z^ 



(30) 



To end this section, recall that we have throughout assumed M > 0. When M < 0, the 
transformation ()2ip becomes ill defined, but it can be verified that the transformation from 
()20p directly to (I24p remains well defined. For M = 0, ()24p is the extremal spinless BTZ 
hole with the angular direction unwrapped, and the Killing horizon of dt is degenerate and 
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has vanishing temperature. For M < 0, (j24p has a massive point particle at r = with the 
angular direction unwrapped, and dt does not have a Killing horizon. It is therefore only with 
M > 0, or in the limiting case M = 0, that we can associate to the metric ()24p a temperature 
and an entropy. 



4 Energy-momentum and thermodynamics 

Our metric ()20p has the asymptotically AdS form of ()15h and ()16p . where the conformal 
boundary metric is the Milne metric (jl7p . On the conformal boundary of AdSs, this Milne 
universe covers the diamond in which the region U < 0, y > of (I22p meets the conformal 
boundary [23]. 

The boundary energy-momentum tensor can be calculated using holographic renormaliza- 
tion equations. In our case of a two-dimensional boundary, the formula is [19] 



T,. = ^y'J - (31) 



Working in the co-moving Milne coordinates (r, rj) of (jlZp , we find 

_ C{M - 1) / r-^ 0\ 

~ IGttGs [ I )■ ^^^^ 

While (j32p should now contain the energy-momentum tensor of the fluid whose thermo- 
dynamics we wish to recover, it could also contain a vacuum energy contribution. As the 
temperature (j27p of the bulk black hole vanishes in the limit M — > 0, we expect this to be the 
limit in which the boundary fluid disappears. We therefore split the total energy-momentum 
tensor as T^^, = T^IJ'^^'^'* -|- T^T'^^ , where 

- ^{\' (33) 

rir' = -TT^i '„ " ]• (34) 

and we interpret T^!!"''^^ and T^^^^^ as respectively the fluid and vacuum contributions. 

As T^^i/ ^ is not proportional to the boundary metric, it cannot be the energy-momentum 
tensor in any Poincare-invariant vacuum state, not even after any Poincare-invariant renor- 
malisation. It is however invariant under the boosts generated by 9^, and it is proportional 
to the energy-momentum tensor of a massless scalar field in the conformal vacuum of the 
Milne universe [24] . The boundary vacuum state from which T^™'^'' arises is therefore uncon- 
ventional from the ion collision perspective but instead adapted to the conformal invariance 
of the gauge/gravity duality. 

The fluid contribution to the energy-momentum tensor takes the perfect fluid form, 

Tit"'^ = {e + p)u^u,+pg^;^l (35) 
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where the fluid's normahsed velocity vector is u'^ = (1,0) and the energy density e and 
pressure p are given by 

The fluid is thus comoving in the Milne universe, following the (1 + l)-dimensional version 
of the Bjorken flow ([1]). 

This was straightforward; the real problem is to define for this time dependent situation 
a time dependent temperature T{t) and entropy density s(r) so that s{T{t)) = dp{T)/dT 
with p defined in (I36p . We are not aware of a method that could be justified with the same 
rigor as in the static case, but we now present two independent arguments, both of which 
lead to the same result. 

• The first argument starts from the entropy density. Recall that in the static case the 
entropy density was given by (j28p with Cdrj as the longitudinal volume element. For 
the time dependent expanding case the longitudinal volume element is, instead, rdi]. 
Thus the proper entropy density and temperature (the latter from s ^ T for an ideal 
gas in 1 + 1 dimensions) are 

Inserting ([371) ([36]) . we find 

IT L TT L 

<r) = p{t) = j^T^r), s{t) = ^^(r). (38) 

Note that this result satisfies the proper thermodynamic relation s{T) =p'[T). 

• The second argument starts from the temperature. The conformal boundary metric in 
()30p is the Minkowski metric, dscBTZ ~ —dt^ + C'^drf , and it follows from (|2ip and 
(|23p that T / C = e*/'^. The boundary metrics in (|2Up and (j3Up are thus related by the 
time-dependent conformal scaling 

f^^Miine = -dr^ + r'dv' = {r/Cf [-df + C^dri") = {r/Cf d^BTZ- (39) 

In a time-independent scaling of a static metric, the temperature scales as inverse 
distance, as seen from the period of thermal Green's functions. If this scaling is assumed 
to hold also in our time-dependent situation, equations (j27p and (j39p lead again to the 
time-dependent temperature and entropy density ()37l) . the entropy density now being 
fixed by the proportionality argument s ~ T. 

Finally, recall that the gravity side of the duality is not just AdSs but AdSaxSs x T4. The 
parameter C/G3 = 4Qi(55 is fixed by ([5]), and we arrive at the final result 

e(T) = p{T) = vrQiQsT^, s{T) = 2ttQiQ^T. (40) 



8 



We may compare this with (1 + l)-dimensional ideal gas in thermal equilibrium with Nb = 
^QiQ5 bosonic and Nf = 4Q1Q5 fermionic massless degrees of freedom, for which 

e{T)=p{T) = {Nb + ^Nf)^T^ = TrQiQ5T^ s{T) = 2t:QiQ^T. (41) 

The results duly coincide. For AdSsxSs the boundary pressure can be computed in the 
strong coupling limit g^N(. ^ 1. There its value is 3/4 times the weakly coupled ideal gas 
value. 

To end this section, we wish to compare our results for the boundary thermodynamics to 
a method that was applied to the corresponding (3 + l)-dimensional situation in [Tj. Let us 
assume M > 1 and compare the time-dependent metric (|20p and the static metric (j30p . The 
latter looks like the former if one replaces Zh = 2£/\/M by = vt, where v is given by (I25p . 
and also replaces C'^drf by r'^drf'. This suggests that the static temperature Tbtz = ^/{T^Zh) 
([27ll should be replaced by the time-dependent temperature 



VM- 1 



(42) 



■KZfi TTVT 2irT 

For the entropy, we can derive a consistent result using the area formula 

' = •* = /*'V^ = /'*)|2r = H£/rd,. (43) 

where 7 is the determinant of the metric on the one-dimensional hypersurface of constant r 
and z. To convert S to an entropy density, we have to divide by the volume, which again is the 
standard longitudinal boost invariant expression / rdr]. The entropy density thus becomes 



C 1 £ VM^l_ (44) 



2G3 VT 4G3 r 
Using formula (j36p for the energy density and the pressure, we find 

.M=pM = ^(l-±)-'r^M, s(r) = ^T(r). (45) 

The proper thermodynamic relation s(T) = p'{T) is thus recovered in the semiclassical limit, 
M ^> 1. Alternatively, if we view the total boundary energy-momentum tensor ()32p as coming 
from a fluid, without a vacuum component, formula (j36p is modified by the replacement 
M ^ M — 1, and the thermodynamic relation s{T) = p'{T) is recovered from ()42p and (j44[l 
for ah M > 1. 

The result from our method for the temperature and entropy density is given by Eq. (j37p 
and holds for M > 0. This result and the method of [7] are thus in agreement in the 
semiclassical limit, M ^> 1. In this limit the temperature satisfles vrTr = \/M/2 ^> 1 and 
is therefore consistent with the uncertainty principle. We view this agreement as indirect 
support for the semiclassical limit of the thermodynamical conclusions obtained in [7], where 
the global structure of the (4-|- l)-dimensional bulk spacetime remained open. We emphasise, 
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however, that our method is coordinate- invariant: the bulk solution (j20p defines a specific 
way of extending the boost isometry on the boundary into a spacelike translational isometry 
in the bulk, and in our method the thermodynamics arises from the Killing horizon of the 
timelike Killing vector that commutes with these spacelike translations. From the geometric 
viewpoint, we therefore view our method as more reliable beyond the semiclassical limit, at 
least in the regime in which a classical bulk solution without quantum corrections can be 
expected to give accurate results for the boundary quantum theory [251 126] . 



5 Bulk metric with rotation 

If the starting ansatz (llOjl is generalised to include a term proportional to drdr], integration of 
the field equations with our boundary condition (jlTj) yields one more constant of integration, 
corresponding to the rotation of the black hole. We shall now find this rotating generalisation 
of (|20p by working backwards from the rotating BTZ metric \17\ I18j. 

ds"^ = -Fdt'^ + ^ + r^{dip + N^dtf, (46) 
F 

where 

This metric satisfies Einstein's equations ([7]) and is therefore locally isometric to AdSa. The 
parameters r± determine the mass parameter M and the angular momentum parameter J 

by 

M = J = (48) 

When If is periodic with period 2tt and the parameters satisfy < |r_ | < r_|_, or equivalently 
< |J| < CM, the metric (|46p is the nonextremal BTZ black hole. The Boyer-Lindquist 
-type coordinates {t,r,(p) are singular at r = r-|_, which is the Killing horizon of the Killing 
vector dt, and if r_ 7^ 0, there is also a coordinate singularity at r = |r_|, which is an inner 
Killing horizon of dt- Both of these horizons are nondenegenerate. For r_ = the metric 
reduces to (j24p . and the conformal diagram was shown in Figure [TJ For r_ 7^ 0, the conformal 
diagram is shown in Figure [2j The ADM energy and the angular momentum at the infinity 
r — > 00 are respectively M/{^G^) and J /{SG^). For us is not periodic, but the crucial point 
is again that the observations about the Killing vector dt remain valid even for nonperiodic ip, 
and the only Killing vectors that commute with translations in </? are c?<^ and dt- We shall 
continue to describe the metric in the black hole terminology even for nonperiodic (p. 

For definiteness, assume for the moment that M > 1 and (M — 1)^ — [J / C)^ > 0, which 
is a subcase of the nondegenerate black hole range \J\ < CM. Starting from the exterior 
r > r^, we first transform from (t,r,ip) to (/?, r, ry) by 

t = C[(3- f{r)], ^ = rj + cj' N^{r)f'{r)dr, (49) 
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Figure 2: The conformal diagram of the non-extremal BTZ black-and-white hole spacetime with 
J ^ [T7l|T8]. The solid (respectively dotted) lines are curves of constant r (constant z) in the region 
< z < VT oi the metric ([ST]), where M > 1, (M — 1) — (J/C) > and v is given by (f53|) . r increases 
upwards and z increases to the left. The dashed line is the coordinate singularity z = vt. 
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where / is a function of r only and satisfies 



The metric becomes 



dr 



+ r^{dr] + CN^dpy 



We then define a as the positive solution of 



{r/cy 



+ e 



+ 



M - 1 



where 



It follows that 



F + l 



(M - 1)^ - 
dr 



{j/cy 



1/4- 



da. 



e-" + e-2« + ii;2(M-l) 



p2a 



(50) 
(51) 

(52) 
(53) 

(54) 
(55) 



Substituting ([52]) . and ([55]) in ([5T]) yields in the coordinates (/?, a,r/) a metric that is 
regular for a > 0, rather than just for r > r+. Finally, defining the coordinates (r, ry, z) by 



z/C 



ve 



,/3-a 



./3 



(56) 



the metric takes the form 



ds' 



(M-l)z2 1 r 



2r2 



+ 



16 



+ 1 + 



(M 



-l)z^ l_ 
2t2 + 16 



(M- 1)^- (J/£)' 
(M - if - {J/Cf 



dr' 



Jz^ 



dr dr] 



-Kr-'drf + dz'^ 



(57) 



This is the promised rotating generalisation of (f20l) . reducing to ((20]) for J = 0. 

Although the transformations shown above assumed M > 1 and (M — 1)^ — {J / Cf > 0, an 
analytic continuation argument (or an explicit computation) shows that ()57p solves Einstein's 
equations ([7]) for any values of M and J. One of the Killing vectors is djj, and the metric 
has at z — > the asymptotically AdS form of (I15p and (I16p with the conformal boundary 
metric ()17p . The transformation to the BTZ form ()46p can be found for any M and J by a 
suitable generalisation of the above formulas. 

When < |J| < CM, so that the spacetime is a nondegenerate black hole, the region 
covered by the metric (I57p can be found as in section [3l Consider in particular the case 
in which M > 1 and (M — 1)^ — {J / Cf > 0. The coordinate singularity in (I57p is then at 
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z = VT, and a regular region that is asymptotically AdS at z ^ is < z < vt, corresponding 
to Q > in the coordinates {p,a,rj). Examination of the coordinate transformations given 
above and the properties of the extended BTZ spacetime \17\ [T8] shows that the situation is 
qualitatively similar to that with J = 0. The coordinate singularity is again on a spacelike 
surface in the white hole region, and this coordinate singularity is neither an event horizon, 
apparent horizon nor a dynamical horizon [211 122j . The curves of constant r and the curves 
of constant z are shown in the conformal diagram in Figure [2l 



6 Energy- momentum from the rotating bulk 



The boundary energy- momentum tensor for the bulk metric ()57p . with arbitrary values of M 
and J, can be computed directly from (j3ip . In the coordinates {T,rj), the boundary metric 
is the Milne metric (1171) and we obtain 



C ( M-l -{J/C)t 



167rG3r2 I -(J/£)r (M - l)r2 



(58) 



As the components in (jSSp do not depend on 77, the energy-momentum tensor is invariant 
under the boosts generated by the boundary Killing vector 9^. This had to be the case since 
dr^ is a Killing vector also in the bulk. 

As in section HI we split the energy- momentum tensor as T^jy = T^^^'^^ + T^^^ , where the 
vacuum contribution T^™^^ (|34p is what remains in the limit of vanishing M and J and the 
fluid contribution T^^"''^'' is 

,„ .,^ 1 / CM -Jt \ 

167rG3r2 1^ -Jt CMt^ ) ' ^ ' 

A nonzero J clearly affects T^^^^^ . We wish to understand how. 

When |J| < £|M|, T^^'^^^ ^ can be written in the perfect fluid form (j35p . where the 
normalised velocity vector of the fluid is 



w 



(cosh (p,T ^ sinh 0) , (60) 



(p is given by 

tanh2(^ = -4T, (61) 
and the energy density e and pressure p of the fluid in its local rest frame are 



CM / J Y 

As the components in (j60p do not depend on rj, the velocity vector field is invariant 
under the boosts generated by 0^, and u'^ has at each point the peculiar velocity Vpec = 
tanh (p relative to the flow field d-^ of the co-moving Milne observers. From the gauge theory 
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viewpoint, the boundary energy-momentum tensor therefore describes a perfect fluid flow 
that is invariant under the longitudinal boosts and is at each point moving with respect to 
the Bjorken similarity flow ([1]) with the longitudinal velocity Upec- Note, however, that the 
flow lines are not inertial when ^ 0- In the Minkowski null coordinates {x~^,x~), in which 
x^ = re^''/\/2 and ds'^ = —2dx~^ dx~ , we have 




and the flow lines are 



.- = e-V^. (63) 



A'(i-)"''<'*', (64) 



where the positive constant K labels the lines. It can be verified that on each line the proper 
time A is given by A = ^/2x~^x~ / cosh(i;^)) and the proper acceleration has the magnitude 
|tanh(^|/A. While all the flow lines start from {x'^^x') = (0,0) at A = 0, the proper 
acceleration diverges as A ^ 0, and although the proper acceleration approaches zero as 
A ^ oo, the fall-off is so slow that the flow lines are not asymptotically inertial at A — > oo. 

If, in addition to \J\ < C\M\, we assume also M > 0, the bulk solution is a black hole, and 
we can use its thermodynamics to equip the flow (j64p with a temperature and an entropy 
density as in section SI It remains however an open problem to identify a microscopic gauge 
theory process whose macroscopic properties the flow (I64p and its associated thermodynam- 
ical quantities would describe. 

For completeness, consider briefly the remaining ranges of the parameters. When \J\ = 
C\M\ 7^ 0, the bulk black hole is extremal and T^l^™'^'' takes the null dust form 

^(fiuid) _ 

where /c^ = (1, ^t^^) is a null vector and the sign in the second component is that of CM/ J. 
When J = CM = 0, T^!^"''^^ vanishes. When CM < \J\, T^!^"'*^^ has no real eigenvectors 
and does not arise from conventional matter fields [21]. We conclude that Tj^^^'^^ is that of 
a perfect fluid with a positive energy density precisely when |J| < £|M|, that is, precisely 
when the bulk solution is a nondegenerate black hole. 



7 Quasinormal modes 

As the metric ()20p with M > does not cover the black hole horizon, but does cover part of 
the white hole region, it is tempting to interpret the ion collision on the boundary as dual to 
an eruption from a white hole in the bulk. We emphasise that our derivation of the boundary 
temperature and entropy density ()37p did not rely on such an interpretation but operated 
directly on the Killing horizon of the extended bulk solution and its thermodynamical prop- 
erties. From a general relativistic viewpoint, one is indeed inclined to read little into the 
behaviour of a specific set of coordinates in the deep bulk region, even when the coordinates 
are near the boundary adapted to the boundary physics of interest: it is always possible to 
introduce a smooth coordinate transformation that is the identity in a neighbourhood of the 
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boundary but drastically changes the region covered deep in the bulk. We show in the Ap- 
pendix that coordinates similar to those in (I20p can be introduced even in (l + l)-dimensional 
Minkowski spacetime, with the coordinate singularity on a spacelike curve in the past light 
cone of the origin. 

All that being said, suppose one did wish to adopt the eruption from a white hole as a 
serious physical picture in the bulk. Would this picture have any observable consequences 
in the boundary physics? We shall now argue that the bulk quasinormal modes \27\ [28l [29] 
may provide insight into this question. 

For concreteness, we consider a bulk scalar field $ that satisfies the Klein-Gordon equation 

O--^)<^ = 0, (66) 



£2 

where ^ > 0. For technical simplicity (cf. [29]), we further assume that yT+T^ is not an 
integer. 

Consider the quadrant U < 0, y > of the spinless BTZ spacetime (I22p in the Schwarzschild- 
like coordinates (p^ . We look for a solution to ([66]) that is independent of r/ (corresponding 
to boost invariance on the boundary) and has the separable form $ = e~*'^*i?(r), where uj is 
a nonvanishing complex number. A pair of linearly independent solutions for R{r) is [27] 



where 



^'iW = f 1 - ^) (^)"'-^("± + 7,^± + 7;2''± + i;i-^), (67) 
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- (i - yr+71) , (69) 



F is the hypergeometric function, and we are assuming that 2i'± + 1 is not a negative integer. 
We write $±(t,r) = e~*'^*i?±(r). Using ([23]) to go to the Kruskal coordinates ([22]) . we see 
that <!>+ continues regularly across the past branch of the horizon at r = r_|_, into the white 
hole region of the spacetime, but it is singular on the future branch of the horizon and cannot 
thus be regularly continued into the black hole region. Conversely, <!>_ continues regularly 
into the black hole region but not into the white hole region. 

We now require ^± to vanish at r ^ oo. Using equation 15.3.6 in [30], and the technical 
assumption that yT+7i is not an integer, we find that this happens precisely when iy± = 
— (n + 1) -|- 7, where n = 0, 1, . . .. Dropping an overall constant, the solutions then take the 
form 



^±,nit,r) = exp 



2r+(n + l-7)t 
£2 



„2 \ / ^2 \ 7-(«+l) 




xFUn, -n; 2(1-7); , n = 0,l,.... (70) 
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^-,n are the usual bulk quasinormal modes [27], decaying in t by falling into the black hole. 
These modes are singular at the white hole horizon, but this singularity in the past does not 
pose a problem for using $-,n to describe decay processes whose initial conditions are set in 
the exterior. By contrast, the time-reversed bulk modes $+,n erupt from the white hole and 
become singular on reaching the black hole horizon. Because of this singularity, $_|_^^ are not 
usually considered relevant for bulk physics whose initial conditions are set in the exterior. It 
seems however less clear whether this singularity would preclude $+,n from describing physics 
that takes place on the boundary, since the black hole horizon is not in the causal past of the 
boundary. 

Thus, if the eruption from a white hole is proposed to have a physical meaning as the dual 
to the boundary ion collision, a possible consequence is that the relevant solutions to the 
wave equation (|66p should be the eruption modes, rather than the usual quasinormal modes. 
From (12ip and (j23p we find that ^>±,n have at z — > the asymptotic form 



The quasinormal modes therefore have on the boundary a decreasing power-law behaviour 
in T, as one expects of a relaxation process. However, the eruption modes have an increasing 
power-law behaviour in r for large n, and even for all n if the hole is so large that r+/£ > 1. 

Finding a power-law increase in some thermodynamical variables on the boundary would 
thus provide smoking-gun evidence for the eruption picture in the bulk, and it would also 
suggest a similar picture in the physically more interesting system with a (3-|- l)-dimensional 
boundary [7J . Conversely, the absence of a power-law increase in the boundary termodynamics 
would discourage the eruption picture as one with physical consequences for the boundary. 
We leave further scrutiny of this question to future work. 

8 Conclusions 

Gauge theory/gravity duality has made interesting predictions about matter in static thermal 
equilibrium [Il|2]. An obvious problem is to explore what, if anything, the same framework 
can say about systems in expansion. 

In this paper we have studied this problem for (1 + l)-dimensional matter that expands 
as the Bjorken similarity flow. Using an ansatz (jlOp adapted to the symmetries we found an 
explicit time dependent AdSa bulk solution (j20p with a time dependent coordinate singularity 
in the bulk. By explicit coordinate transformations we showed that the solution is part of the 
spinless BTZ black hole, with the angular dimension unwrapped. The coordinate singularity 
is on a spacelike surface in the white hole region. 

The holographic energy-momentum tensor on the boundary was computed with standard 
techniques and separated into two components, a vacuum contribution (coming from the 
bulk metric with M = 0) and a contribution that corresponds to the boundary fluid. This 
separation turned out important for a consistent definition of the temperature and entropy 
density of the fluid. The fluid was found to be an ideal gas in adiabatic expansion. We also 
recovered the known fact that the boundary pressure, as calculated in the AdSsxSs x T4 




(71) 
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approach, is exactly the same as that of a massless ideal gas of an appropriate number of 
bosons and fermions in one spatial dimension. 

To obtain the time-dependent temperature and entropy density of the fluid, we first com- 
puted the time-independent temperature and entropy density of the BTZ hole, defined in a 
standard way with respect to the longitudinal volume element £dr/, and we then performed 
the time-dependent scaling £ ^ r to the longitudinal volume element rdr] that is appropriate 
for the boost symmetry of the Bjorken flow on the boundary. We emphasise that the thermo- 
dynamic results therefore relied in no way on the coordinate singularity in the metric (j20p . 
In a more geometric language, the bulk solution ()20p defines a specific way of extending the 
boost isometry on the boundary into a spacelike translational isometry in the bulk, and the 
thermodynamics arose from the Killing horizon of the timelike Killing vector that commutes 
with these spacelike translations. 

Our AdSa solution, of course, does not contain the full dynamical content of the physically 
interesting case where the boundary is (3 + l)-dimensional. For example, in our case e(r) is 
proportional to for all r, while on the (3+l)-dimensional boundary e(r) appears to encode 
qualitatively different physics at large [7] and small |14j r. However, in our case the global 
structure of the bulk metric is completely known. We were in particular able to verify that 
our coordinate-invariant derivation of the thermodynamics was in the semiclassical regime 
(M ^ 1) fully compatible with the method of [7], in which the temperature and entropy 
formulas of a static black hole are extended to the time-dependent case by comparing the 
coordinate singularities. 

We also showed that inclusion of angular momentum in the bulk leads to a boundary 
flow that is still boost invariant but has a nonzero longitudinal velocity with respect to the 
Bjorken expansion. Finally, we argued that the bulk quasinormal modes may shed light on 
the possible physical relevance, or lack thereof, of the observation that the coordinates in 
(j20p cover part of the bulk white hole interior but none of the black hole interior. 

The boundary flow temperature (j37p can be written in the form 

T{T)=Ti^, ^Tin = \VM. (72) 
r z 

Interesting physical questions here are what is the smallest time Tj for which (j72p holds, what 
is the thermalisation time, and what is the associated maximum temperature. To address 
these questions, additional information about the process would need to be introduced. This 
information could regulate the singularity in ([72]) by 1/r 1/(t + tq) or produce other 
corrections, for example of the type 1/r — > l/r — tq/t"^ . As the gauge theory /gravity duality 
approach is designed to give vacuum expectation values of gauge theory operators, it is by 
no means obvious how it should be used for processes starting at some r = 0. 

Acknowledgements. We thank Esko Keski-Vakkuri, D. Podolsky, Eric Poisson and Jung- 
Tay Yee for discussions. This research has been supported by Academy of Finland, contract 
number 109720 and by PPARC grant PP/D507358/1. JL thanks Helsinki Institute of Physics 
for hospitality during the final stages of the work. 



17 



Appendix: Singular coordinates in (1 + 1) Minkowski 



In this appendix we give an example of a coordinate system in (1 + l)-dimensional Minkowski 
spacetime with a singularity structure similar to that in (j20p . 

Let {T,X) be the usual Minkowski coordinates, ds'^ = —dT"^ + dX"^ . We introduce the 
coordinates (r, z) by 



T-X 
T + X 



\d.{vt/z) + 1 



\[x].{vt/ z) — 1] r, 

where < z < oo, < r < oo, and t> is a positive constant. The metric becomes 



^.2 = 4 



^)ln(!^ 



dT^ + dz^ 



(73) 



(74) 



This metric is regular for Q < z < vt and < ?;t < z, and the oordinate singularity at 
z = VT is on the spacelike curve T = -V^2 + 1. The curves of constant r and the curves of 
constant z in the region Q < z < vt are qualitatively similar to those shown in Figure [TJ 
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